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The magneto-oscillations in graphene bilayers are studied in the vicinity of the K and K' points 
of the Brillouin zone within the four-band continuum model based on the simplest tight-binding ap- 
proximation involving only the nearest neighbor interactions. The model is employed to construct 
Landau plots for a variety of carrier concentrations and bias strengths between the graphene planes. 
The quantum-mechanical and quasiclassical approaches are compared. We found that the quantum 
magneto-oscillations are only asymptotically periodic and reach the frequencies predicted quasi- 
classically for high indices of Landau levels. In unbiased bilayers the phase of oscillations is equal 
to the phase of massive fermions. Anomalous behavior of oscillation phases was found in biased 
bilayers with broken inversion symmetry. The oscillation frequencies again tend to quasiclassically 
predicted ones, which are the same for K and K', but the quantum approach yields the gate-tunable 
corrections to oscillation phases, which differ in sign for K and K' . These valley-dependent phase 
corrections give rise, instead of a single quasiclassical series of oscillations, to two series with the 
same frequency but shifted in phase. 

PACS numbers: 71.20.-b, 71.70.Di, 81.05.ue 



I. INTRODUCTION 

In solids subject to a magnetic field B, the energy spec- 
trum of charge carriers is quantized into Landau levels 
(LLs). The magneto-oscillations (MOs) observed in the 
Shubnikov-de Haas and de Haas-van Alphen effects re- 
flect the oscillations of the density of states (DOS) with 
the field intensity. The DOS reaches maxima at magnetic 
fields, B n , for which the LLs with the index, n, cross the 
Fermi energy, Ep. 

The Landau plot is a plot of inverse magnetic fields, 
1/B n , versus the LL index, n. It is a standard tool used 
to determine the frequency and phase of MOs, and the 
related important characteristics of the investigated sys- 
tems. 

The construction of the Landau plot is based on the 
Onsager-Lifshitz quasiclassical quantization rule^^ 
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where A(Ep) is an area of the extremal cross-section of 
the Fermi surface (FS) cut by the plane perpendicular 
to the magnetic field direction, e is the electron charge 
and 7 is a constant which describes the phase of MOs. 
It follows from Eq. ([1]) that MOs of DOS are periodic in 
1/B and their frequency F is related to A(Ep) by 
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The Onsager-Lifshitz quantization rule has been origi- 
nally designed for three dimensional metals, where the 
validity of the quasiclassical approximation is guaran- 
teed by a large number of LLs bellow Ep in accessi- 
ble magnetic fields. However, the method is also widely 
used when two-dimensional (2D) systems are investi- 
gated. Here, the importance of F is stressed out by the 



fact that the carrier concentration is proportional to the 
area surrounded by the Fermi contour. 

In general, the rule should not be applicable to 2D 
systems. Subject to strong magnetic fields, the quan- 
tum limit with only one LL below Ep can be easily 
reached. But in the majority of such systems, the pe- 
riodicity of MOs is preserved due to the simple parabolic 
(Schrodinger-like) energy spectra of the 2D electron lay- 
ers in the semiconductor structures, which yields the LL 
energies proportional to B. 

In 2004 a single sheet of graphene was separated from 
bulk graphite by micromechanical cleavage 4 It was con- 
firmed experimentally that electrons in graphene obey a 
linear energy dependence on the wave- vector k, as pre- 
dicted many years ago by the band structure calculation. 4 
Both electron and hole charge carriers behave like mass- 
less relativistic particles - Dirac fermions (DFs), and 
there is no gap between the valence and conduction 
bands. The electron and hole Dirac cones touch at a 
neutrality point. 

Subject to a magnetic field B, the DFs form LLs with 
energies proportional to \fB. In the seminal papers^^ 
the Shubnikov-de Haas MOs in graphene were found pe- 
riodic in 1/B, similarly as in the 2D gas of Schrodinger 
fermions (SFs) with the parabolic energy spectra, but 
with the phase shifted by ir. The shift, which was clearly 
demonstrated by the Landau plot of magneto-resistance 
oscillations, is due to the existence of the zero-energy LL 
in the linear Dirac spectrum, shared by electrons and 
holes. Note that 7 = 1/2 for SFs, and 7 = for DFs. 

In addition to a single layer graphene, also a few layer 
graphene samples can be prepared. Among them a bi- 
layer graphene (BLG), in which two carbon layers are 
placed on top of each other with a standard Bcrnal stack- 
ing, is of particular interest. Probably the most re- 
markable feature of this structure is the possibility to 
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open a gap between the valence and conduction bands 
through the application of an external field or by chemi- 
cal doping. 7-9 . This phenomenon is closely related to the 
gate-induced breaking of the inversion symmetry of the 
crystal J£r— 

Note also that the application of the gate voltage is a 
necessary condition for the experimental observation of 
MOs in BLG. Without a gate voltage, the sample is neu- 
tral, the Fermi energy is located in the neutrality point, 
and no free charge carriers should be present in perfect 
samples. 

There are two ways of how to apply the gate voltage. If 
the external voltage is applied symmetrically from both 
sides of a sample, just Ep and the concentration of car- 
riers are varied, and no gap is opened. The tunable gap 
appears in the presence of external electric field resulting 
from the asymmetrically applied gate voltage. 

Let us point out that the charge carriers in BLG are 
neither SFs nor DFs, and therefore it is of interest to 
construct the corresponding Landau plots to see how far 
the bilayer energy spectra from these two simplest possi- 
bilities are. 

This task is simplified by the fact that the electrochem- 
ical potential (i.e., also Ep) is kept constant during mag- 
netic field sweeps in gated samples. According to Ref. LLJ 
and [TH carrier density oscillations are compensated by 
gate current oscillations in the case of fixed Ep. Note 
that in bulk samples, where the charge neutrality must 
be preserved, the carrier concentration is considered to 
be fixed. 

To construct the Landau plot, we will first calculate 
the quasiclassical frequencies of MOs in BLG, based on 
their zero-magnetic-field electronic structure. 

Later on we will compare these quasiclassical frequen- 
cies with results of the quantum-mechanical calculation 
of the electronic structure of BLG subject to a perpen- 
dicular magnetic field. 



II. ZERO-FIELD ELECTRONIC STRUCTURE 

The electronic structure of BLG can be described by 
the simple tight-binding model involving only the nearest 
neighbor interactions 

A single layer honeycomb lattice, with two atoms per 
unit cell, results from two superimposed triangular lat- 
tices labeled A and B. The unit cell is defined by the 
lattice vectors a\ and 0,2, making the angle 60°, the lat- 
tice constant a is equal to 2.46 A. The bilayer is formed 
by two graphene sheets, 1 and 2, arranged in the Bernal 
stacking. The distance between layers is 3.37 A. Thus the 
unit cell of a bilayer has four atoms, its lattice structure 
is sketched in Fig. [TJ 

In addition to the intralayer parameter 70 and the in- 
terlayer parameter t, the corresponding Hamiltonian de- 
pends on the potential energy difference between the two 
layers, which we denote 2u. The parameter 70 3.1 eV 
yields the Fermi velocity vp « 1.0 x 10 6 m/s, defined 




FIG. 1. (Color online) Lattice structure of a graphene bilayer. 
The unit cell is a green parallelepiped. 



by hvp = jQ^/3a/2. We further consider that t « 0.39 
eV, and the energy 2u varies between and 250 meV. 23 
While 70 and t are fixed by nature, we assume that u and 
Ep are the adjustable parameters. 

If we employ the continuum approximation, 4 the 
Hamiltonian H in the vicinity of the K point can be 
written as 



I Hb 2 b 2 Hb 2 a 2 Hb 2 a 1 Hb 2 b 1 \ 
Ha 2 b 2 Ha 2 A 2 Ha 2 Ax Ha 2 b x 



H = 
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where the matrix elements of the first layer are given by 



Ha 1 a 1 — Hb 1 b 1 
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Similarly, the matrix elements corresponding to the sec- 
ond layer read 



Ha 2 a 2 — Hb 2 b 2 — u, 

Ha 2 B 2 = = Kv F {k x 



3 2 A 2 — !<>"*■ \nx -r 

There are only two nonzero interlayer matrix elements 
H Al A 2 — Ha 2 a 1 = t. 

The Hamiltonian H' in the vicinity of the K' point 
has a similar structure, the matrix elements of H' are 
complex conjugates of the matrix elements of H. 

The above Hamiltonians can be diagonalized 
analytically. 16 ' 17 i 20 ' 22 i 24 The zero-field energy branches 
of the conduction band, E c i(k) and E C 2(k), and the 
valence band, E v \{k) and E V 2(k), of a bilayer result 
from hybridization of Fermi cones of layers 1 and 2, 
mediated by the interlayer matrix element t. Note that 
E v i(k) = —E c i(k) and E v2 (k) = —E c2 (k) and that the 
valley degeneracy is preserved, i.e., we get the same 
bands in valleys K and K' . 
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Fermi contours 



FIG. 2. (Color online) The „mexican hat" shape of the va- 
lence and conduction bands of a biased bilayer. The blue and 
red colors correspond to higher probability of finding charge 
carriers in the layers 1 and 2, respectively. Three groups of 
the Fermi contour are possible depending on the value of Ep: 
the double circles (Al), the circles (A2), and the Fermi rings 
(A3). 



For u = two Fermi cones are replaced by four bond- 
ing and antibonding hyperbolic bands. The bonding va- 
lence and conduction bands, E v \(k) and E c i(k), touch 
at k = 0, the separation between bands of a bonding- 
antibonding pair is equal to t on the energy scale. 

When the interlayer voltage is applied, the Fermi cones 
of two layers are shifted along the energy axis, and the 
separation of the neutrality points becomes equal to 2u. 
The hybridization due to the interlayer parameter t is 
strongest near the cone cross-points. The resulting four 
bands are shown in Fig. [2] It turns out that for any 
finite u a gap is open between the topmost valence band 
E v i(k) and the bottom conduction band E c i(k). The 
conduction band acquires a ,,mexican hat" shape with 
energy minima at nonzero k and a local maximum at 
k = 0. We can write 



(0) = u, 

(0) = V~i 



rp max 
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E c f n (k) = A = ut/y/4v? + t 2 . 



t 2 , 
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Note that for large k the band E c x(k) describes electrons 
localized mostly in the layer 1. Near the local maximum 
at k = the holes in the layer 2 prevail. Similar conclu- 
sions can be drawn for the topmost valence band E v \(k). 

As mentioned in Introduction, the quasiclassical fre- 
quencies of the bilayer, F\ and F2, are proportional to 
areas surrounded by the Fermi circles, which depend, for 
a given u, on the Fermi energy value Ep. Three different 
possibilities are depicted in Fig. [2] for the case of conduc- 
tion/valence bands. (For the valence bands Ep should 
be replaced by —Ep.) 

The analytic expressions for the quasiclassical frequen- 
cies -Fi and F 2 read 
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the frequencies Fi(2) are even functions of variables Ep 
and u. The frequency F2 is equal to zero at the local 
maximum E t ^ iax (Q), and at the minimum i? c ™ m (0). For 
a finite u, the frequency F2 approaches F\ at E^ m (k). 

Three forms of the Fermi contour are possible depend- 
ing on the value of Ep. First, the large Ep cuts both 
conduction bands and Fi > F2 > 0. The frequency 
F\ corresponds to electron orbits localized mainly in the 
layer 1, the frequency F 2 corresponds to hole orbits lo- 
calized mainly in the layer 2. Second, only the band E c \ 
is cut by Ep. Then F\ > and F 2 < 0. In this case F 2 
is just a parameter and does not have the meaning of a 
true frequency. At last, the Ep cuts the bottom conduc- 
tion band E c i(k) twice, if it is less than a local energy 
maximum, i? c ™ ax (0). Then again Fi > F 2 > 0. In that 
case F\ is the frequency of an electron orbit in the layer 
1 while F2 belong to a hole orbit in the layer 2. Close 
to the local minima the difference between electrons and 
holes is smeared and charge carriers are present in both 
layers as indicated by the change of line colors in Fig. [5] 

For the special case of u = 0, Eq. ((5]) reduces to 
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Then the gap between the valence and conduction bands 
as well as the local maximum i? c ™ ax (0) all disappear. 

The quasiclassical phases of MOs are not accesible via 
the Onsager-Lifshitz quantization rule, Eqs. ([T]) and ([2]). 
To find the energy spectra beyond the quasiclassical ap- 
proximation, we need to diagonalize the magnetic Hamil- 
tonians H and H'. 



III. MAGNETIC FIELD EFFECTS 

The magnetic Hamiltonians can be obtained from the 
zero-field Hamiltonians by modification of matrix ele- 
ments H AlBl , H Bl Ai , H A2 b 2 and H B2 a 2 The ma- 
trix elements of the magnetic Hamiltonian in the vicinity 
of the K point are 



H AlBl =H BiAi = ^2\e\hv F Bn, 

H A2B2 = H* B2A2 = ^2\e\hv F B{n + \). 

The other matrix elements remain the same as in the 
zero-field Hamiltonian. Near the K' point, 



H 'a 1 b 1 



2\e\hv 2 F B(n + l), 



A 2 B- 



- H B 2 A 2 



2\e\hv 2 F Bn. 



We need not diagonalize these Hamiltonians to con- 
struct the Landau plot. If we look for magnetic fields 
B n at which the LLs cross Ep, it is enough to find the 
poles of the resolvent G(z) = (z — as it defines the 

density of states g(Ep) through 



(5) 



1 



g(E F ) cx lmTrG(E F + i0). 



(7) 
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The easiest way to find the poles is to solve the corre- 
sponding secular equation for B n assuming the fixed Ep. 

We start with the simplest case of the unbiased BLG 
(u = 0). Then the secular equations can be given a very 
convenient form, utilizing the quasiclassical frequencies 
of MOs, presented in the previous paragraph, Eq. (|6]), 

B 2 n{n + l)-B(n+^\ (Fl + F 2 ) + FiF 2 = 0. (8) 

The Hamiltonians H and H' yield identical equations for 
valleys K and K'. 

While the secular polynomial is quartic in energy it is 
only quadratic in B. Therefore, to construct the Landau 
plot it is enough to solve the quadratic equation to find 
B n in terms of fixed E = Ep,. 

The quasiclassical phase 7 can be easily obtained from 
Eq. ([5]). For a large number of LLs below Ep one may 
assume that n(n + 1) — > (n+ 1/2) 2 , and then Eq. ((5J can 
be written in the form 



B n 



1 



-Bin 



1 



(Fi + F 2 ) + FiF 2 = 0. (9) 



From here we obtain the asymptotic quasiclassical Lan- 
dau plots 



(2) 
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i.e., we found that the phases of MOs correspond to SFs 
with 7 = 1/2, in agreement with quasiclassical treat- 
ments of systems with inversion symmetry. Note that F 2 
is positive only in the rather unrealistic case \Ep\ > t. 

To get Landau plots for an arbitrary n we can express 
the solution of Eq. (IS]) as 



2F X F 2 1 1 // IV 4F!F 2 
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or, if we define dimensionless 5 by 



S = 



( F\ — F 2 



\F 1 +F 2 J ' 
we can write (see also Ref. [27 
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Here the negative sign in the numerator/denominator 
corresponds to the frequency Fj in the quasiclassical 
limit, and the positive sign to the quasiclassical frequency 
F 2 . It is obvious that for 8 7^ the MOs are not periodic 
in 1/B. 

The case of the biased BLG (u ^ 0) must be treated 
separately, as the presence of the electric field perpendic- 
ular to layer planes breaks the inversion symmetry and 
lifts the valley degenerac y 11 ' 12 



The secular equations can be given again a form 
quadratic in B, but the coefficients do not depend ex- 
clusively on the quasiclassical frequencies as in Eq. ((8|). 
We can write 



B 2 n(n+ 1) - B 



{n+^)(F 1 +F 2 ) + F 



F 1 F 2 = 
(14) 

for the Hamiltonian H in the vicinity of K. In compari- 
son with Eq. (|8]) there is an extra term 
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In the vicinity of K' we obtain a very similar equation 
from the Hamiltonian H ' , the only difference is that Fo is 
replaced by — Fo. The extra term, Fo, is the reason of the 
valley asymmetry. It is obvious that Eq. (fT4)) gives two 
different series of solutions, B n , for positive and negative 
F . 

The quasiclassical frequencies Fl and F 2 are even func- 
tions of Ep and u. It means that there are the same fre- 
quencies not only for K and K', but also for the electrons 
and holes with energies Ep and — Ep, respectively. Note 
also that Fi and F 2 do not depend on the sign of u. On 
the other hand, Fq is an odd function of Ep and u. Thus 
F) breaks the K - K' symmetry, and also the symmetry 
between the electron and hole oscillations with the same 
quasiclassical frequencies. The change of sign of u also 
reverts the roles of K and K' valleys, i.e., what is valid 
for K with u > is valid for K' with u < 0. 

Also Eq. (|14p can be rewritten to an equation similar to 
Eq. (|13[) , but with an additional dimensionless parameter 
A, which depends on Fq, 



A 



"0 



Fi +F 2 

Then the analytic solution reads 



(16) 



(F+F 2 ) 2 

( U ) F m n+\ + \T^{n+\ + \y + n(n + 1)(1 - 5) 



B n 



1 ± 



(17) 

This equation reduces to Eq. (fT3"|) for A = 0. 

Now it is a more difficult task to find an asymptotic 
expression for the oscillation phase than in the previous 
case u = 0. If we solve Eq. (fT4|) for n + 1/2 we get 



1 F + F 2 ± yj (Fi - F 2 ) 2 + 4FF + B 2 
n+ 2= 2B ' (18) 



which for B approaching zero yields 

F U2) , 1 c 



(19) 



Here £ is a gate-tunable correction to the oscillation 
phase, given by 



F, 



Epu 



Fi-F 2 y/ (E F — u 2 )t 2 + AEpu 2 ' 



(20) 
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This correction differs in sign for K and K' and also 
differs for electrons and holes from the same valley with 
the same absolute value of energy. 



IV. RESULTS AND DISCUSSION 

In the unbiased BLG the energy u is equal to zero 
and the parameter 6, which appears in Eq. (1131) . has a 
particularly simple form 



F 



(21) 



For small Fermi energies only the bottom branch 
E c i (k) of the conduction subband is cut by Ep and only 
the frequency F\ is defined. For Ep approaching zero, 
the parameter 5 diverges. This implies that for energies 
close to the band bottom Eq. (fT3"]) can be written as 
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the form found for the extremal electron and hole orbits 
in graphite^! which clearly indicates the aperiodicity of 
oscillations. 

The Fermi energies greater than t are rather unrealis- 
tic. Nevertheless we can consider this hypothetical case 
in our theoretical treatment. We can write, for Ep = t 
and 6 = 1, Eq. (|13p as follows 



El 
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The Landau plots calculated for two selected values of 
Ep, Ep < t and Ep > t, which cross the dispersion 
curves are presented in Fig. [3] The Landau plots are 
the same for Ep and — Ep due to the inversion sym- 
metry conservation. One can observe that in the unbi- 
ased bilayer the phases of MOs, corresponding to both 
frequencies F\ and F 2 , approach the phase of massive 
fermions,7 = 1/2, for higher quantum numbers of LLs. 

In BLG, an applied electric field leads to asymmetry 
between K and K' valleys that gives rise to nontrivial os- 
cillation phenomena in magnetic fields. To illuminate an 
anomalous behavior of oscillations, we plotted in Fig. [4] 
the field dependence of LLs in BLG. 

In a single layer graphene the LL fans of electrons and 
holes start at the zero-field neutrality point. The neu- 
trality points of two independent layers are shifted by 2u 
and the LLs of holes from the layer 1 cross the LLs of 
electrons from the layer 2, as shown in Fig. [4] by thin 
brown lines. 

In BLG the shape of LL spectrum results from hy- 
bridization of LL spectra of layers 1 and 2. Due to the 
interlayer interaction, represented by the matrix element 
t, we have four fans of LLs which start at zero-field en- 
ergies E v2 (0), E vl (0), E cl (0) and E c2 (0). 
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FIG. 3. (Color online) (a) The electronic bands of the unbi- 
ased graphene bilayer (w = 0). The horizontal lines denote 
the Fermi energies which cross the electron and hole disper- 
sion curves, (b) The ,,phases" 7i(2)(^f) = F\m\/B — n, for 
Ep depicted in a), plotted as functions of the Landau level 
index n. 



FIG. 4. (Color online) The electron and hole Landau levels 
(in the K' valley) of two layers are mixed by the interlayer 
interaction, t. For the energy range corresponding to Fermi 
rings in zero magnetic filed (see A3 in Fig. [2]), Ef cuts the 
Landau levels twice. This is the reason for the anomalous 
phase in the quasiclassical limit B — > 0. 
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u = OmeV(K K) 
u-O.OSeV(K) 
u = 0.05eV(K') 
u = 0.125 eV (K) 
u = 0.125 eV (K-) 



FIG. 5. (Color online) The „phase" 71 = Fj/B — n calculated 
for the fixed quasiclassical frequency F\ = 70 T and various 
u, for the electron K and K' valleys as a function of the LL 
index, n. 



The hole levels from the layer 1 and the electron levels 
from the layer 2 avoid to cross, and the low-field hole LLs 
smoothly turn into the electron LLs as B increases. This 
is indicated in Fig. 2] by the change of LL color from 
red to blue. The LLs from a fan starting at zero-field 
energy E c \ (0) have minima in their field dependence and, 
therefore, can be cut twice by a single Ep. Moreover, the 
minima are not the same for all levels and, consequently, 
not all levels are cut by a single Ep. 

This is reflected in the quasiclassical approach as the 
gate-dependent correction to the MO phase, £, which is 
related to the energy difference 2u between two layers. 
Note that in the region of energies corresponding to the 
Fermi rings the expression ([20)) diverges at E c 1 lln (k) and 
is equal to 1/2 for Ep = E™ a *(0). As the many body 
effects can play a role in this low concentration range, the 
above one-electron picture is probably oversimplified. 
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FIG. 6. (Color online) The DOS of the unbiased (a) and 
biased (b, c) BLG versus dimensionless value of the Lan- 
dau plot, Fi/B, with the fixed quasiclassical frequency F\ = 
70 T. The frequency F\ corresponds to the situation when 
F\ > 0, F2 < 0, i.e., only the lowest conduction/valence en- 
ergy band is cut by Ef- In (b, c) the blue peaks show DOS 
calculated for the K valley, whereas the red ones are related 
to the K' valley. 



We start our discussion with the simplest case, 
E c i/vi(Q) < \Ef\ < E c2 / V 2{0), when only one conduc- 
tion/valence energy band is cut by Ep (see Fig. [2]). The 
single quasiclassical frequency F% corresponds to a single 
Fermi area, which is the same for K and K' . 

According to Eq. (|19p . the electron peaks in DOS are 
shifted by £ to the left in the K valley, whereas the peaks 
in the K 1 valley are shifted by £ to the right. The shift 
magnitude ranges from zero to 1/2 depending on the en- 
ergy difference between two layers. In FigHJthe ,, phase" 
71 is plotted as the function of n for three different cases 
with u equal to 0, 0.05 and 0.125 eV. The Fermi energies 
are chosen to keep the same Fermi area (and the same 
fixed F±) in all three systems. Only for u = the curves 
are identical for K and K', for it ^ the curves are 
substantially different. 

In Fig.[6]the shifts of peaks in the above three cases are 
shown explicitly. There is a single series of oscillations for 
the unbiased bilayer, as the valley degeneracy is preserved 
in a system with the inversion symmetry. 

The effect of the gate-tunable valley splitting originates 
in two series of oscillations which differ for the different 
u. Let us emphasize that all series of oscillations have the 
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FIG. 7. (Color online) (a) The electron/hole bands of the 
biased graphene bilayer with the gap 2u — 0.25 eV at k = 0. 
The horizontal lines denote the Fermi energies which cross 
the electron (solid lines) and hole (dashed lines) dispersion 
curves, (b) The ,, phases" 71(2) (Ef) = Fuas/B — n calculated 
for Ef depicted in (a) plotted as functions of the LL index, 
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FIG. 8. (Color online) The same as in Fig. only for the 
lowest electron/hole bands and A < \Ef\ < u. 



same quasiclassical frequency F\, but the quasiclassical 
phases depend on the choice of u and on the choice of 
the valley index. 

We complete our discussion with cases when the Fermi 
energy cuts the conduction/ valence bands twice. The 
Landau plots of the biased bilayer with u = 0.125 eV, 
which is probably the highest experimentally accessible 
valued calculated for four selected Fermi energies, two 
in the conduction band, and two in the valence band, 
are presented in Fig. [7] In accordance with types of the 
Fermi contours in Fig. [2j the first and second cases are 
depicted. 

The situation is more complicated in the region of 
energies, A < \Ep\ < u, for which Ep cuts the lowest 



subband E c \/ V \(k) twice, which is characteristic for the 
third type of the Fermi contours, as shown in Fig. [5] 
The bottom of E c \[k) is at ~ 0.105 eV. The parameter 
£ is far from the values expected for the phase of quasi- 
classical oscillations, it reduces/grows heavily when Ep 
approaches the bottom of E c \(k) / E v \(k) . For Ep — u, 
£ becomes closer to —1/2 for Ep in the conduction band 
and 1/2 for Ep in the valence band. 

V. CONCLUSIONS 

Using a four-band continuum model, we calculated an- 
alytically the Landau plots in biased and unbiased BLG 
subject to external perpendicular magnetic fields. 

It turns out that the magneto-oscillations are only 
asymptotically periodic, and that in the unbiased bilay- 
ers their phase is equal to the phase of massive fermions. 
The convergence to the quasiclassical limit is slow, and 
depends strongly on the the value of Ep. The conver- 
gence is slower for higher values of Ep. 

Anomalous behavior of oscillation phases was found 
in biased bilayers with broken inversion symmetry. The 
oscillation frequencies again tend to quasiclassically pre- 
dicted ones, which are the same for K and K', but the 
quantum approach yields the gate-tunable corrections to 
oscillation phases, which differ in sign for K and K' . 
These valley-dependent phase corrections give rise, in- 
stead of a single quasiclassical series of oscillations, to 
two series with the same frequency but shifted in phase. 

We also found that for Ep in the region of energies 
corresponding to the Fermi rings in the quasiclassical ap- 
proach, only a limited number of LLs can cut the Fermi 
energy and thus a limited number of magneto-oscillations 
can be achieved. Moreover, their quasiclassical phases 
reach very large values. As the many body effects can 
play a role in the corresponding concentration range, the 
above one-electron picture is probably oversimplified. 
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